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Abstract 

We review some of the most recent results on the dynamics of the Hamiltonian 
Mean Field (HMF) model, a systems of N planar spins with ferromagnetic infinite- 
range interactions. We show, in particular, how some of the dynamical anomalies 
of the model can be interpreted and characterized in terms of the weak-ergodicity 
breaking proposed in frameworks of glassy systems. We also discuss the connections 
with the nonextensive thermodynamics proposed by Tsallis. 
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1 Introduction 



Metastability, nonextensivity and glassy dynamics are features so ubiquitous 
in complex systems that are often used to characterize them or, more in gen- 
eral, to define complexity [1] . In this work we will discuss such features in the 
context of the so called Hamiltonian Mean Field (HMF) model, a system of 
inertial spins with long-range interaction. The model, originally introduced 
by Antoni and Ruffo [2], has been thoroughly investigated and generalized in 
the last years for its anomalous dynamical behavior [3,4,5,6,7,8,9,10,11,12,13]. 
With respect to systems with short-range interactions, the dynamics and the 
thermodynamics of many-body systems of particles interacting with long- 
range forces, as the HMF, are particularly rich and interesting. In fact, more 
and more frequently nowadays, the out-of-equilibrium dynamics of systems 
with long-range interactions or long-term correlations has shown physical sit- 
uations which can be badly described within the ergodic assumption that is 
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at the basis of the Bohzmann-Gibbs thermostatistics. hi all sTich cases it hap- 
pens, for instance, that a system of particles kept at constant total energy E, 
does not not visit all the a-priori available phase space (the surface of con- 
stant energy E) , but it seems to remain trapped in a restricted portion of that 
space, giving rise to anomalous distributions that differs from those expected. 
A few years ago, Tsallis has introduced a generalized thermodynamics for- 
malism based on a nonextensive definition of entropy [14]. This nonextensive 
thermodynamics is very useful in describing all those situations characterized 
by long-range correlations or fractal structures in phase space [6,15,16]. On the 
other hand, the latter feature is also connected with the so called "weak ergod- 
icity breaking" scenario, which is at the basis of the long-term relaxation and 
aging observed in glassy systems. Such systems show competing interactions 
(frustration) and are characterized by a complex landscape and a hierarchical 
topology in some high dimensional configuration space [17], which, in turn, 
generates a strong increase of relaxation times together with metastable states 
and weak chaos. 

The Hamiltonian Mean Field model, considered in this paper, is exactly solv- 
able at equilibrium and exhibits a series of anomalies in the dynamics, as the 
presence of quasistationary states (QSS) characterized by: anomalous diffu- 
sion, vanishing Lyapunov exponents, non-gaussian velocity distributions, ag- 
ing and fractal-like phase space structure. Furthermore, the model is easily 
accessible by means of both molecular dynamics and Monte Carlo simulations. 
Thus, it represents a very useful "laboratory" for exploring metastability and 
glassy dynamics in systems with long-range interactions. The model can be 
considered as a minimal and pedagogical model for a large class of complex 
systems, among which one can surely include self-gravitating systems [13] and 
glassy systems [10], but also systems apparently belonging to different fields 
as biology or sociology. In fact, we recently found similar features also in the 
context of the Kuramoto Model [19], one of the simplest models for synchro- 
nization in biological systems [20]. Moreover, the prohferation of metastable 
states in the vicinity of a critical critical point in the phase diagram seems to 
be responsible for the onset of complexity and diverging time calculation in 
many different kind of algorithms [18]. 

In this paper we focus on two different aspects of the HMF model: its glassy- 
dynamics and the possible connections with the generalized thermodynamics. 
The paper is divided into two parts. In Section 2.1 we investigate the model 
following the analogy with glassy systems and the 'weak ergodicity breaking' 
scenario. In previous works we have shown that the "thermal explosion" , char- 
acteristic of initial conditions with finite magnetization, drives the system into 
a metastable glassy- like regime which exhibits 'dynamical frustration'. With 
the aim to characterize in a quantitative way this behavior, we have explic- 
itly suggested to introduce a new order parameter, the 'polarization', able to 
measure the degree of freezing of the rotators (or particles). Here we present 
new numerical results reinforcing the glassy nature of the QSS's metastability 
and the hierarchical organization of phase space. In Section 2.2 we investigate 
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the links with nonextensive thermostatistics. In rcf. [11] wc have found that, 
for a particular class of initial conditions with constant velocity distribution 
and finite magnetization, the velocity correlations obtained by integrating the 
equations of motion of the HMF model are well reproduced by q-exponential 
curves. Here, we show numerical evidences that the superdiffusion observed 
in the anomalous QSS regime (ref. [5]) can be linked with the q-cxponential 
long-term decay of the velocity correlations, as analitically suggested by a for- 
mula obtained by Tsallis and Bukman [21] for a non/mear Fokker-Planck (FP) 
equation, using an ansatz based on the generalized entropy. 



2 Anomalous dynamics in the HMF model 

The HMF model has been introduced originally in ref. [2] with the aim of 
studying clustering phenomena in A^-body systems in one dimension. The 
Hamiltonian of the ferromagnetic HMF model is: 

N 2 1 JV 

H^K + v^Y.^ + —j:[i-cos{e,-ej)] , (1) 

i=l i,i=l 

where the potential energy is rescaled by 1/A^ in order to get a finite spe- 
cific energy in the thermodynamic limit — > oo. This model can be seen as 
classical Xy-spins (inertial rotators) with unitary masses and infinite range 
coupling, but it also represents particles moving on the unit circle. In the lat- 
ter interpretation the coordinate 6i of particle i is its position on the circle 
and Pi its conjugate momentum (or velocity). Associating to each particle the 
spin vector 

^ {cos 6i, sin 6i) , 

it is possible to introduce the following mean-field order parameter: 

M^^IE^.I , (2) 
-'^ i=i 

representing the modulus of the total magnetization. 

The equilibrium thermodynamical solution in the canonical ensemble predicts 
a second-order phase transition from a low-energy condensed (ferromagnetic) 
phase with magnetization M 7^ 0, to a high-energy one (paramagnetic), where 
the spins are homogeneously oriented on the unit circle and M = 0. The caloric 
curve, i.e. the dependence of the energy density U = H/N on the temperature 
T, is given hy U = |- + |(1 — M'^) [2,5]. The critical point is at energy density 
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[/c = f , corresponding to a critical temperature Tc = |. 

At variance with the equihbrium scenario, the out-of-cquiUbrium dynamics 
shows, just below the phase transition, several anomalies before complete equi- 
libration. More precisely, if we adopt the so-called Ml initial conditions (i.e.), 
i.e. 9i — for all i (M(0) = 1) and velocities uniformly distributed {water bag), 
the results of the simulations, in a special region of energy values (in particular 
for 0.68 < U < Uc) show a disagreement with the canonical prediction for a 
transient regime whose length depends on the system size N. In such a regime, 
the system remains trapped in metastable states (QSS) with vanishing mag- 
netization at a temperature lower then the canonical equilibrium one, until 
it slowly relaxes towards Boltzmann-Gibbs (BG) equilibrium, showing strong 
memory effects, correlations and aging. This transient QSS regime becomes 
stable if one takes the infinite size limit before the infinite time limit. 



2. 1 Dynamical frustration and hierarchical structure 

As required by the discovery of correlations and aging, which in turn imply 
complex trajectories of the system in phase space, it is interesting to explore 
directly the microscopic evolution of the QSS. This can be easily done plotting 
the time evolution of the Boltzmann fi-space, where each particle of the system 
is represented by a point in a plane characterized by the conjugate variables 
Oiandpi, respectively the angular position and the velocity of the ith particle . 
It has been shown[8] that, during the QSS regime, correlations, structures 
and clusters formation in the /x-spacc appear for the Ml i.e., but not for ini- 
tial conditions with zero magnetization, the so called MO i.e.: in the latter 
case both the angles and velocities distributions remain homogeneous from 
the beginning and a very slow mixing of the rotators has been observed. For 
the Ml case, the dynamics in //-space can be clarified through the concept 
of "dynamical frustration": the clusters appearing and disappearing on the 
unit circle compete one with each other in trapping more and more particles, 
thus generating a dynamically frustrated situation that put the system in a 
glassy-like regime. 

In Fig. 2.1 we show a molecular dynamics simulation where the complete distri- 
bution function f{9,p,t) is considered for different values of time. In fact, we 
plot - for Ml i.e., N=10000 and U=0.69 - a sequence of snapshots of f{6,p,t) 
for six different times: at the beginning of the simulation (t=0), in the QSS 
regime (t=50-500) and towards canonical equilibrium (t= 10000). In the QSS 
region one clearly observes the presence of competing clusters, each cluster be- 
ing composed by particles with both angles and velocity included in the same 
/i-space cell (notice that, in our simulations, we considered a total of 100x100 
cells for the /i-space lattice). For t=10000, instead, any trace of macroscopic 
glassy behavior has disappeared. 

In Fig. 2 we show the power law behavior of the cluster size cumulative dis- 
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Fig. 1. Snapshots of f{9,p,t) for t=0 (upper left), 50 (upper right), 100 (center 
left), 200 (center right), 500 (lower left) and 1000000 (lower right). In this case we 
considered N=10000 at U=0.69 and Ml i.e.. See text. 

tributions calculated in the case U = 0.69 for several snapshots in the QSS 
regime at time t=200,350 and 500. For each one of the 100x100 cells a sum 
over 20 different realizations (events) has been performed. Then, for each clus- 
ter size (greater than 5 particles) the sum of all the clusters with that size has 
been calculated and plotted. As one can see from 2, the distribution does not 
change significativcly in the plateaux region as expected. We report also a 
power law fit (drawn as a straight dashed line above the data points) which 
indicates that the cluster distribution has an approximately exponent decay 
— 1.6. The cluster size distribution reminds closely that of percolation at the 
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Fig. 2. We plot for the case U = 0.69 and N = 10000 (using Ml i.e.) the cumulative 
distribution of the clusters size, calculated for three times {t = 200, 350, 500) in the 
QSS region. See text. 

critical point, where a lenght scale, or time scale, diverges leaving the system 
in a self-similar state [22]. More in general, it has been also suggested [23] that, 
optimizing Tsallis' entropy with natural constraints in a regime of long-range 
correlations, it is possible to derive a power-law hierarchical cluster size dis- 
tribution which can be considered as paradigmatic of physical systems where 
multiscale interactions and geometric (fractal) properties play a key role in the 
relaxation behavior of the system. Therefore, we can say that the power-law 
scaling resulting in the distributions of Fig. 2 strongly suggests a non-ergodic 
topology of a region of phase space in which the system remains trapped dur- 
ing the QSS regime (for the Ml i.e.), thus supporting the weak-ergodicity 
breaking scenario. A weak breakdown of ergodicity, as originally proposed by 
Bouchaud for glassy systems [24] , occurs when the phase-space is a-priori not 
broken into mutually inaccessible regions in which local equilibrium may be 
achieved, as in the true ergodicity-breaking case, but nevertheless the system 
can remain trapped for very long times in some regions of the complex energy 
landscape. In fact it is widely accepted that the energy landscape of a finite 
disordered (or frustrated) system is extremely rough, with many local min- 
ima corresponding to metastable configurations. Since the energy landscape is 
rough, these local minima are surrounded by rather high energy barriers and 
we thus expect that these states would act as "traps" which get hold of the 
system during a certain time r. 

In ref . [24] such a mechanism has been proposed in order to explain the aging 
phenomenon, i.e. the dependence of the relaxation time on the history of the 
system, i.e. on its age t^- Actually, it results that r^ax — tw, being r^ax the 
longest trapping time actually encountered during a waiting time t^- In other 



6 



words, the deepest state encountered traps the system during a time which 
is comparable to the overall waiting time, a result that - in turn - allows to 
quantitatively describe the relaxation laws observed in glassy systems [17]. 
Aging phenomenon has been found also in the HMF model for Ml i.e., more 
precisely in the autocorrelation functions decay for both the angles and ve- 
locities [9] and for velocities only [8], thus reinforcing the hypothesis that a 
weak ergodicity-breaking could really occur in the metastable QSS regime and 
could be related to the complex dynamics generated by the vanishing of the 
largest Lyapunov exponent and by the dynamical frustration due to the many 
different small clusters observed in this regime. Such a scenario is in agree- 
ment also with the results about anomalous diffusion shown in ref. [4]. where 
the probability distribution of the trapping times, calculated for a test particle 
in the transient QSS regime for Ml i.e., shows a clear power law decay 

Ptrap ~ t~'' , (3) 

characterized by an exponent u related to the anomalous diffusion coefficient. 
In the next section we will show that the anomalous diffusion coefficient can 
in turn be connected with the velocity correlations decay by means of the 
nonextensive formalism, thus suggesting a deeper link between the latter and 
the weak ergodicity breaking. 

2.2 Nonextensive thermodynamics and HMF model 

In previous works it was shown that the majority of the dynamical anomalies 
of the QSS regime, among which /i-space correlations, clusters and dynamical 
frustration, are present not only for Ml initial conditions, but also when the 
initial magnetization M{t = 0) is in the range (0, 1] [11]. 
In order to prepare the initial magnetization in the range < M < 1, we 
distribute uniformly the particles into a variable portion of the unitary circle. 
In this way we fix the initial potential energy V{6) and, in turn, the magne- 
tization. Finally, we assign the remaining part of the total energy as kinetic 
energy by using a water bag uniform distribution for the velocities. The ve- 
locity correlations can be calculated by using the following autocorrelation 
function [11] 

CW = ^Ep.(%(0) , (4) 

where Pj{t) is the velocity of the j-th particle at the time t. In Fig.2.2-left, 
we plot the velocity autocorrelation function (4) for N = 1000 and M(0) ~ 
1,0.8,0.6,0.4,0.2,0. An ensemble average over 500 different realizations was 
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Fig. 3. (Left) Correlation functions vs time for different initial magnetization (sym- 
bols). The solid lines are normalized q-exponentials. (Right) We plot the mean 
square displacement of the angular motion o"^ oc t"' vs time for different initial mag- 
netizations. The exponent 7 which fit the data and charcterize the behavior in the 
QSS regime is also reported for the different initial conditions used. The dashed 
lines have a slope corresponding to these values. 

performed. For M(0) > 0.4 the correlation functions are very similar, while 
the decay is faster for M(0) = 0.2 and even more for M(0) = 0. If we fit these 
relaxation functions by means of the Tsallis' q-exponential function 

e,{z) = [l + {l-q)z]T^ , (5) 

with z = and where r is a characteristic time, we can quantitatively dis- 
criminate between the different initial conditions. In fact we get a q-exponential 
with q = 1.5 for M(0) > 0.4, while we get q = 1.2 and q = 1.1 for M(0) = 0.2 
and for M(0) = respectively. Notice that for q = 1 one recovers the usual ex- 
ponential decay [15,5,6,8]. Thus for M(0) > correlations exhibit a long-range 
nature and a slow power-law decay. This decay is very similar for M(0) > 0.4, 
but diminishes progressively below M(0) = 0.4 to become almost exponential 
for M(0) = 0. 

In order to study diffusion, one can consider the mean square displacement of 
phases cr^(t) defined as 

1 ^ 

^"(t) = -^T.m - mr =< m - > ■ (6) 

where the symbol < ... > represents the average over all the rotators. 
The quantity o"^(t) typically scales as cT^(t) ~ f^. The diffusion is normal 
when 7 = 1 (corresponding to the Einstein's law for Brownian motion) and 
ballistic for 7 = 2 (corresponding to free particles). For different values of 7 the 
diffusion is anomalous, in particular for 1 < 7 < 2 one has superdiffusion. We 
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Fig. 4. In this figure we check the correctness of the 7 — g conjecture (see text) for 
the HMF model by plotting 7 as a function of the ratio 2/(3 — g) for different initial 
conditions with variable magnetization. See text. 



notice that the quantity cr^(t) can be rewritten by using the velocity correlation 
function C{t) as 



a\t) = Jdti Jdt2 < Vjih) Vjiti) > ='^J^^^J Citi) , (7) 
00 00 



where C{t) is defined as in Eq.4. 

Superdiffusion has been already observed in the HMF model for Ml initial 
conditions [4]. Rercently we have also checked that, even decreasing the initial 
magnetization, the system continues to show superdiffusion [11]. We illustrate 
this behavior in Fig.2.2-right, where one sees taht, after an initial ballistic 
regime (7 = 2) proper of the initial fast relaxation, the system shows su- 
perdiffusion in correspondence of the QSS plateau region and afterwards. The 
exponent goes progressively from 7 = 1.4 — 1.5 for 0.4 < M(0) < 1 to 7 = 1.2 
for MO. In the latter case, we have checked that, by increasing the size of the 
system, diffusion tends to be normal (7 ~ 1 for N=10000). The slow decay and 
the superdiffusive behavior illustrated in Fig. 2. 2 can be connected by means 
of a conjecture based on a theoretical result found in ref. [21] by Tsallis and 
Bukman. In fact in that paper the authors show, on general grounds, that non- 
extensive thermostatistics constitutes a theoretical framework within which 
the unification of normal and correlated (driven) anomalous diffusions can be 
achieved. They obtain, for a generic linear force F{x), the physically relevant 
exact (space, time)-dependent solutions of a generalized Fokker-Planck (FP) 
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equation 



by means of an ansatz based on the Tsallis entropy. 

For our purpose, we remind here that such a FP equation, in the nonhnear 
"norm conservation" case (z/ 7^ 1 and ^ = 1), generates Tsalhs space-time 
distributions with the entropic index q being related to the parameter v by 
q — 2 — u. By means of the latter, and following again ref.[21], it is possible to 
recover the following relation between the exponent 7 of anomalous diffusion 
(being oc t"^) and the entropic index q 

2 2 

7 = = (9) 

' 1 + 1/ 3-g ^ ' 



Hence, being in diffusive processes the space-time distributions linked to the 
respective velocity correlations by the relation (7), one could think to investi- 
gate if the relation (9) would be satisfied choosing the entropic index q, char- 
acterizing the correlation decay, and the corresponding anomalous diffusion 
exponent. This is done in Fig.4 where in order to check the latter hypothesis, 
that we call the 7 — ? conjecture, we report the ratio 2/{3-q) exponent 
7 for various initial conditions ranging from M(0)=1 to M(0)=0 and different 
sizes at U = 0.69. Both q and 7 have been taken from the results shown in 
Fig. 2. 2. Within an uncertainty of ±0.1, the data show that this ratio is always 
one, thus providing a strong indication in favor of this conjecture, which is 
satisfied for the HMF model. 

Summarizing, we have shown numerical simulations which connect the su- 
perdiffusion observed in the anomalous QSS regime of the HMF model to the 
g-exponential long-term decay of the velocity correlations in the same regime. 
This new result is very interesting because opens a way to set a rigorous an- 
alytical link between the entropic index q and the dynamical properties of 
nonextensive Hamiltonian many-body systems. 



Conclusions 



We have briefiy reviewed some of the anomalous features oberved in the dy- 
namics of the HMF model, a kind of minimal model for the study of complex 

behavior in systems with long-range interactions. We have also discussed how 
the anomalous behavior can be interpreted within the nonextensive thermo- 
statistics introduced by Tsallis, and in the framework of the theory glassy 
systems. The two pictures are not in contradiction and probably have more 
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strict links than previously thought, which deserve to be further explored in 
the future. 
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